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Single-particle spectra and two-particle correlation functions measured by the 
NA49 collaboration in central Pb+Pb collisions at 158 GeV/nucleon are presented. 
These measurements are used to study the kinetic and chemical freeze-out condi- 
tions in heavy ion collisions. We conclude that large baryon stopping, high baryon 
density and strong transverse radial flow are achieved in central Pb+Pb collisions 
at the SPS. 



1 Introduction 

The goal of studying high energy heavy-ion collisions is to understand the 
properties of nuclear matter under extreme conditions of high energy density 
(> 1 GeV/fm 3 ). It is unlikely under such conditions that spatially separated 
hadrons exist. Instead, formation of quark-gluon-plasma (QGP), in wbkjh 
quarks and gluons are deconfined over an extended volume, is expected. Era€rQ 
This statejpf matter might have existed in the early universe shortly after the 
Big Bang.U 

If a QGP is produced in a heavy ion collision, the collision system will 
undergo transitions from hadrons (initial nucleons in the incoming nuclei), 
to partons, to interacting hadrons, to, finally, freeze-out particles where the 
measurements are realized. Many signatures of QGP formation have been 
proposed, u These include^nhanced production of strangeness, 13 charmH and 
lcptons, El jet quenching, E3 collective radial flow, liil etc. However, none of 
the signatures is conclusive, and many can be produced by secondary parti- 
cle interactions. E3li3 Therefore, systematic studies of multiple observables are 
necessary in the searching of the QGP formation in heavy ion collisions. 

CERN SPS experiment NA49 has measured many observables of Pb+Pb 
collisions at_L58 GeV/nucleon. These measurements include transverse energy 
production, li-3 baryon stopping and negative hadron distributions, E-3 negative 
hadronlla and proton til two-particle correlation functions, and charged and 
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functions are presented in Sect. 
Conclusions are drawn in Sect. 



neutral strange particle distributionsHllin central collisions, and directed and 
elliptic flow in non-central collisions. cEl In this talk, only selected results from 
central Pb+Pb collisions are presented. Distributions of net baryons, negative 
hadrons, and charged kaons are presented in Sect. || to address the chemical 
freeze-out conditions. Negative hadron and proton two-particle correlation 

3] to address the kinetic freeze-out conditions. 
The results which have been published or 
have been submitted for publication are not duplicated here but are referred 
to. All the plots shown here are preliminary. 

NA49 is a large acceptance hadron experiment, measuring about 800 
charged particles in a central Pb+Pb event. The experimental apparatus con- 
sists of four large time-projection chambers (TPCs). Two of them are placed 
along the beam axis inside two dipole magnets, which have maximum inte- 
grated field strength of 9 Tesla-meters. The other two are downstream of the 
magnets on either side of the beam axis. Further downstream are four time-of- 
flight (TOF) walls covering smaller acceptance than the TPCs. The particle 
identification was performed by measuring, the specific ionization (dE/dx) de=j 
posited by charged particles in the TPCs,E3 and the velocities of the particlesca 
in the TOF walls. The dE/dx resolution was 5%, and the TOF resolution was 
60 ps. A beam of Pb nuclei struck on a Pb target of thickness 224 uwJ-cm 2 
placed in front of the first TPC. A zero-degree hadronic calorimeter, tilEJ lo- 
cated downstream along the beam-axis, measured mainly the energy remnant 
in the projectile spectators. By applying a threshold on the measured energy 
as an experimental trigger, the 5% most central events were recorded. These 
events corresponded to collisions with impact parameter below 3.3 fm. 



2 Single-particle distributions 

Extensive studies of Si+A and Au+Au collisions at the BNL AGS and S+A 
collisions at the CERN SPS have shown that nucleus-nucleus collisions are not 
mere superpositions of nucleon-nucleon collisions (N+N) at the corresponding 
energies. U Protons are shiftcdJjux^p: into the central-rapidity region (baryon 
stopping) in these collisions BBEZEa than in N+NEa. The NA49 result on 
the net baryon rapidity distribution in central Pb+Pb collisions over a wide 
rapidity range (Fig. 1 in_R£fll3) shows a larger average rapidity shift than that 
in S+S (3% centrality) .E3c3 This is consistent with the picture of more multiple 
collisions suffered by each nucleon on average in Pb+Pb than in S+S. 

A larger fraction of energy is converted from the longitudinal direction into 
the transverse direction in Pb+Pb central collisions than in S+S. Two conse- 
quences can be expected: more particle production and/or stronger transverse 
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expansion in Pb+Pb than S+S. Products from high energy heavy ion col- 
lisions are dominantly pions. The rapidity distribution of negative hadrons 
(h~), mainly 7r _ 's, scales with number of participants from S+S to Pb+Pb 
central collisions almost perfectly (Fig. 2 in RefE3), in contrast to the different 
net baryon distributions. Since the total pion multiplicity provides a good 
approximation to the total entropy produced, the scaling implies that no ex- 
tra entropy per participant was generated in Pb+Pb central collisions with 
respect to S+S. However, one should note that pion production does not scale 
from N+N to S+S. Ej Collective transverse expansion is reflected in particle 
transverse distributions; heavy particles pick up more transverse momentum 
from the expansion than light particles, exhibiting larger inverse slopes of the 
exponential transverse mass (toj_ = \Jp\_ + to 2 ) spectra. In Fig. |], the inverse 
slopes are-plotted against particle masses (to) foij— thxee systems: Pb+Pb 113 
and S+S B at the SPS, and Au+Au at the AGS.BEl The inverse slope sys- 
tematically increases with particle mass in all three systems, consistent with 
the picture of collective transverse radial flow. The inverse slope is larger in 
heavy system Pb+Pb than light system S+S at similar energies, and is larger 
at high energy (SPS) than low energy (AGS) in systems with similar sizes. 

~~> — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — 
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Figure 1: Mass systematics of inverse slopes of particle m± spectra in the central rapidity 
region for three central collision systems: Pb+Pb (5%) at ^/s = 17.3 GeV, S+S (3%) at 
= 19.4 GeV, and Au+Au (4%) at -/s = 4.74 GeV. The straight line is to guide the eye. 

Kaons, carrying most of the strangeness produced in heavy ion collisions 
at SPS energy and lower, have been extensively studied. NA49 has measured 
charged kaon production at mid-rapidity by the TOF walls and over wide 
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rapidity range by the TPCs. The kaon m± spectra, can be well described by 
single exponential functions. At mid-rapidity, the inverse slopes are measured 
to be 240 ± 4 MeV for K+ and 223 ± 3 MeV for K~, and the yields per unit 
rapidity are estimated to be dN/dy = 29.8 ± 0.9 (K+) and 16.2 ± 0.5 (K~). 

The K + 1 K~ ratio is sensitive to the baryon density. In Fig. ||, the ratio at 
mid-rapidity is plotted against the mid-rapidity radjiixof net-baryon yield over 
total- multiplicity ibr_Si-|- A OlJcj and Au+AuE§E30 central collisions at the 
AGS and S+S StZEa and Pb+Pb@ central collisions at the SPS. The total 
multiplicity at mid-rapidity is approximated by the sum of the net baryon yield 
and 3(7r) at mid-rapidity, where (tt) = ((tt + ) + (tt~))/2 for the AGS data and 
(tt) = (h~) for the SPS data. Note that both the abscissa and the ordinate 
are experimentally measured quantities. However, the motivation to plot this 
way is that the ratio of net-baryon yield over total multiplicity is proportional 
to the average baryon density at freeze-out, assuming that the volume of the 
system scales with the total multiplicity. The K + / K~ ratio in e + e~ at LEPE3 
is also plotted. As seen from the plot, the K + / K~ ratio has strong dependence 
on the baryon density. Note that the baryon density in heavy ion collisions is 
higher at low energies and in large systems. 



AGS Si+AI: E802/859 

AGS Si+Au: E802/859 

AGSAu+Au: E866 

SPS S+S: NA35; NA44 Preliminary 

SPS Pb+Pb: NA49 Preliminary 

LEP eV: Particle Data Book 
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Figure 2: Mid-rapidity K + / K~ ratio as function of the mid-rapidity ratio of net-baryon yield 
over total multiplicity in heavy ion collisions at the AGS and the SPS. Rapidity integrated 
K + /K~ ratio in e+e - at LEP is also shown. The straight line is to guide the eye. 



In Fig. |[ thcjtotal^icld K + / K ratios in heavy ion collisions are compared 



to those in N+NBBeI y at various bombarding energies. The ratio strongly 
depends on the bombarding energy in both heavy ion and N+N collisions. 
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However, the ratios are higher in heavy ion collisions than in N+N at the 
same energies. The dependence in N+N can be understood by the higher 
production threshold for K~ than K + . The higher ratio in heavy ion collisions 
is qualitatively consistent with the following conjecture: each nucleon suffers 
multiple interactions in heavy ion collisions, where subsequent interactions 
happen at lower energy than the primary interaction, producing more K + 'a 
relative to K~'s. The ratio of K + / K~ in heavy ion collisions over that in 
N+N (referred as K + /K~ enhancement) is shown in the insert of Fig. [| 
The enhancement is higher in heavy system than in light system, which is 
consistent with the multiple collision picture. However, the enhancement seems 
independent of bombarding energy in systems with similar sizes. 
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Figure 3: Total yield /K~ ratio (except for Pb+Pb, the mid-rapidity ratio is used) as 
function of bombarding energy for N+N (open points) and heavy ion collisions (filled points). 
The ratios of K+ / K~ in heavy ion collisions over these in N+N are plotted in the insert. 



3 Two-particle correlation functions 

Two-particle correlation functions provide space-time information of the parti- 
cles at freeze-out, Q'CJ which cannot be obtained from single-particle distribu- 
tions. Due to quantum interference, the probability to find two identical pions 
(bosons) with low momentum difference is enhanced. The width of the 2it cor- 
relation function after correction for Coulomb interaction is inversely related 
to the effective size of the pion source. Because of the correlation between 
space-time and momentum of the particles produced in heavy ion collisions 
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at freeze-out, due to the collision dynamics and collective behavior such as 
flow, particles within a limited momentum range are emitted from a limited 
space-time region of the source. The two-particle correlation function therefore 
measures, at small momentum difference, a limited region of the source. 

NA49 has measured 2h~ correlation functions, F.p .Based on a model in- 
corporating both longitudinal and transverse flow,c3c3'c3 the transverse and 
longitudinal radii were determined to be 6t=8 fm at mid-rapidity, dropping by 
2-3 fm to forward rapidity (Fig. 4 in Ref. This size, compared to the 1- 
dimensional Gaussian radius erpb ~ 3 fm of the initial nuclei, implies that 
the collision system has gone through significant expansion before freeze-out. 
The xadius parameters decrease with pair transverse momentum p^ (Fig. 6 in 
Rcf. E3). This behavior was also seen in S+A collisions ejS, and is consis- 
tent with transverse radial flow. Combined with the single-particle transverse 
distributions, the freeze-out temperature (assumed to be the same for all par- 
ticle species) and the transverse radial flow velocity at the source surface were 
estimated to be T = 120 ± 12 MeV and fix = 0.55 ± 0.12, respectively. 

Similar p± dependence of the correlation length has been found in Z° 
hadronic decays in e + e~ annihilation at LEP. By including all final particles, 
the LUND model calculation reproduces the dependence. EH The conjecture is 
that contribution from resonance decays, dominant at low px, enlarges the 
correlation length. The p± dependence of the correlation length in e + e~ may 
also be due to space-time-momentum correlation arising from large formation 
time of high p± particles. Study of the effects of resonance decays and particle 
formation time on the source size measurement aLheavy ion collisions, besides 
the effect of transverse radial flow, is underway. Ej 

The baryon density .plays an important role in the dynamical evolution of 
heavy ion collisions. r^ri To measure baryon density, one needs information on 
the space-time extent of the baryon source, which is not necessarily the same 
as that of the pion's. The space-time extent of the baryon source at freeze- 
out can be measured by the two-proton correlation function. The two-proton 
correlatioiuiunction has very different characteristics from the 2n correlation 
function. Due to Fermi statistics, strong interaction between two protons 
is primarily due to 1/4 S- wave attractive interaction of the spin singlet and 
3/4 P-wave repulsive interaction of the spin triplet. The combined strong and 
Coulomb interactions generate a peak in the two-proton correlation function 
at roughly q mv = ^/—q^q^/2 w 20 MeV/c, where q mv is the momentum mag- 
nitude of either proton in the rest frame of the pair, and q^ is the difference 
of the proton four-momentum vectors. The peak height above 1 is inversely 
proportional to the effective volume of the proton source. 

While 2tt correlation function measurements are common, measurements of 
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two-proton correlation functions are rare. NA49 has measured the two-proton 
correlation function at mid-rapidity. This is the first such measurement in 
Pb+Pb collisions at the SPS. The protons are mostly primordial nucleons, 
since baryonpaJitibaryon pair production contribution is at the level of 10% at 
SPS energy. EjlII The measured two-proton correlation function is shown in 
Fig. |. The peak at q inv w 20 MeV/c is evident and the height is 1.14 ± 0.04. 
The correlation function does not go to zero at qi nv = because of the con- 
tamination of protons from weak decays; there, the effect of the experimental 
momentum resolution is negligible. The deviation of the correlation function 
from 1 at large qi im is currently under investigation. 




q irw (MeV/c) 



Figure 4: Two-proton correlation functions (not corrected for protons from weak decays) 
at mid-rapidity measured by NA49 (points), calculated for a Gaussian source (dotted his- 
togram), and calculated for protons generated by RQMD (solid histogram) and by Venus 
(dashed histogram). 

In order to interpret the proton freeze-out conditions, we employed the 
approach of comparing the measured two-proton correlation function with 
theoretical calculations. Given the proton phase space distribution, the two- 
protpn correlation function was calculated using the Koonin-Pratt Formal- 
ism. E30 The experimental momentum resolution was included in the calcu- 
lation. For the proton freeze-out distribution, we used two approaches. (I) 
Gaussian sources were generated with widths a x = a y = a z = a for the space 
coordinates and at = and a for the time duration. The momentum was 
generated by a thermal distribution of temperature T = 120 MeV. The Gaus- 
sian sources do not have space-time-momentum correlation. (II) Protons were 
generated by two dynamical models, the Relativistic Quantum Molecular Dy- 
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namics (RQMD) model (v2.3)SEJ'EE3 and the Venus model (v4.12).B The 
protons have correlation between space-time and momentum intrinsic to the 
dynamical evolution of the models. The strength of the calculated correlation 
functions for (I) were reduced by 30% to account for the contamination of 
protons from weak decays. In the calculation for (II), the experimental accep- 
tance was applied to the protons, and experimental effects in the accidental 
reconstruction of weak decay products as primary protons were incorporated. 
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Figure 5: Peak height above 1 of measured two-proton correlation functions as function of 
beam energy. Contaminations of protons from weak decays are not corrected. The correction 
would bring up the Pb+Pb point by approximately 30%. The corrections are negligible for 
data at lower energies. 



Comparison between the calculations and the experimental data showed 
that the Gaussian source with a x = a y = a z = a t — 3.8 fm well describes 
the data. The calculated correlation function is plotted in Fig. [|. We note 
that a = 3.8 fm is slightly larger than the initial 1-dimensional nuclei size 
dpb rs 3 fm, and is significantly smaller than the radius extracted from the 2h~ 
correlation function. The 1-dimensional effective size of the Gaussian source, 
calculated from the mean separation between protons in the pair rest frame, is 
<xa/v2 = 4.0 fm. We have also studied Gaussian sources with extreme shapes: 
a x = Gy >> a z and a x — a y « a z (both a t = 0). They cannot describe 
the data satisfactorily. The calculated two-proton correlation functions for 
RQMD and Venus are also shown in Fig. |[ Venus overpredicts the strength of 
the correlation function, while RQMD slightly underpredicts its strength. The 
mean separation between protons in the pair rest frame are (ctax, fAy, caz) = 
(6.76,5.34,8.35) fm for RQMD, and (4.41,4.50,5.44) fm for Venus, where z 
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is the longitudinal coordinate. The 1-dimensional effective size, cta/v^ = 
y&Ax&AyCAz/ v2, are 4.74 fm and 3.37 fm, respectively for RQMD and Venus. 
We note that no simple relation exists between a a / \/2 and the apparent single- 
proton source size in the models due to the dynamics. 

We now comment on our two-proton correlation function measurement in 
the context of other measurements. The_nk>p-iSOurce volume is measured to 
be proportional to the pion multiplicity, ol±lj which increases steadily with 
beam energy in similar colliding systems. Lj Due to the large pion-nucleon 
cross-section, one would expect that protons and pions freeze-out under similar 
conditions, therefore, the proton source size would increase with beam energy 
as well. Howevej^.our measurement, inxonjunction with preliminary results 
obtained at GSlBEl and AGS energies, shows that the peak heights of the 
correlation functions in Au+Au and Pb+Pb collisions are rather insensitive 
to beam energy (Fig. [5]). This means that the effective volumes of freeze-out 
protons are similar in these collisions over the wide energy range. More detailed 
studies are needed to understand the instrumental effects and the change in 
space-time-momentum correlation as function of the beam energy. 

4 Conclusions 

We conclude the following for central Pb+Pb collisions at the SPS: 

• Strong baryon stopping. Baryon stopping in Pb+Pb central collisions 
is stronger than in S+S. The baryon spatial density achieved in Pb+Pb 
central collisions is higher than in S+S at the SPS, and is lower than 
in Au+Au at the AGS. The large baryon density resulted in a large 
K + / K~ ratio. The enhancement in the ratio in heavy ion collisions of 
similar systems over N+N does not depend on the bombarding energy. 

• Strong transverse radial flow. Inverse slopes of m<r distributions in the 
central rapidity region are larger in Pb+Pb central collisions than in 
S+S at the SPS and in Au+Au at the AGS. The radius parameters 
from 2h~ correlation function suggest strong expansion. The values of 
freeze-out temperature and transverse radial flow velocity are, however, 
model-dependent. 

• Secondary particle rescattering in heavy ion collisions may be the driving 
force for the large baryon stopping, the strong transverse radial flow, and 
the K + /K~ enhancement. 

• No excessive entropy production. The negative hadron rapidity distri- 
bution scales with number of participants from S+S to Pb+Pb central 
collisions. No indication of excessive entropy production in Pb+Pb cen- 
tral collisions, compared to S+S, is found. 
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• Freeze-out source size. Effective size of the pion source extracted from 
the measured 2h~ correlation function is 6-8 fm at mid-rapidity. Effec- 
tive size of the proton source extracted from the measured two-proton 
correlation function at mid-rapidity is 4 fm. 
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